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Abstract. In the context of the star formation through the fragmentation of an ex- 
tremely metal-deficient protogalactic cloud, the gravitational collapse of filamentary 
gas clouds is explored with H2 and HD chemistry. It is found by ID hydrodynamical 
simulations that the cloud evolution is prescribed mainly by the initial density 
(no) and H2 abundance (a?H 2 ,o)- In particular, it turns out that the evolution of 
low-density filaments (no < 10 5 cm -3 ) bifurcates at a critical H2 abundance of 
' £H 2 ,cr ~ 3 x 1CP 3 , beyond which HD cooling overwhelms H2 cooling. The numerical 

results indicate that the stellar IMF is likely to be double-peaked and deficient in sub- 
solar mass stars, where the high mass peak of the IMF is around WMq or 10 2 Mq, 
dependently on the initial density and H2 abundance. If the gas in protogalactic 
clouds is photoionized by UV radiation or shock-heated, the H2 abundance could 
exceed a;H 2 ,cr — 3 x 10 -3 by H~ reactions. Then, the high mass peak would be 
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Star formation in the early universe should have proceeded in very 
metal-deficient gas. When the metallicity is lower than Z <J 10~ 2 Z & , 
I the thermal property of such metal-deficient gas is essentially the same 

O !~ as that of the metal-free gas because cooling by heavy metals is less ef- 

fective than primordial molecular cooling. Among the primordial molecules, 
the importance of H2 on the star formation in the early universe have 
been emphasized by many authors because they are the most abundant 
molecules in metal-free gas. H2 molecules provide significant cooling 
through rotational and vibrational transitions which can lower the gas 
temperature down to a few hundred K and accordingly reduce the Jeans 
masses to stellar mass scales. 

HD molecules are considered to be the second most abundant molecules 
in primordial gas during the post-recombination era. In spite of its 
low abundance ([HD/H2] ~ 10 -3 ), HD can provide significant cooling 
at low temperature gas because HD has higher radiative transition 
probabilities and lower excitation temperatures than H2. Therefore, 
HD cooling can lower the gas temperature down to T ^ 100K and 
accordingly, the Jeans mass is also reduced. 

In this paper, we examine the effects of the HD cooling on the 
formation of stars in very metal-deficient gas, and elucidate the role 
of HD molecules for the stellar initial mass function (IMF) there. 
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2. Numerical Model and Methods 

Our numerical model and method are the same as those of Nakamura 
and Umemura (2001a). We consider an infinitely long cylindrical gas 
cloud. We deal with the following 14 species: e, H, H + , H~, H2, H^, 
He, He+, He++, D, D+, D~, HD, and HD+. 

The density is assumed to be uniform along the cylinder axis and 
the radial distribution is expressed as p = po (l + r 2 /Rq) , where 

i?o = y/^fkTo/^Gpop) is the effective radius, po is the central mass 
density, To is the initial gas temperature, p is the mean molecular 
weight, and / is the ratio of the gravitational force to the pressure 
force. Our model is specified by four parameters: no, To, /, and the 
H2 abundance £h 2 ,o- It should be noted that a higher xh 2 ,o can be 
translated to a higher electron abundance of a parent cloud. This is 
because when the ionization degree of the parent cloud is higher, the 
H2 molecules form more efficiently via radiative reactions of H + e — > 
H - + hv and H + H~ -> H 2 + e. 

We take into account the following thermal processes: (1) H cooling 
by radiative recombination, collisional ionization, and collisional exci- 
tation, (2) H2 line cooling, (3) cooling by H2 collisional dissociation, 
(4) heating by H2 formation, and (5) HD line cooling. 

3. Numerical Results 

In this section, we examine the collapse of the filaments, taking into 
account the H2 and HD cooling. As shown below, there is a critical 
initial H2 abundance, above which HD cooling predominantly regulates 
the cloud evolution. Also, it is found that the HD cooling does not 
play an important role for high-density gas. Thus, the evolution of the 
primordial filaments is classified into three cases, depending upon the 
initial density and initial H2 abundance; (1) low-density filaments with 
high H2 abundance, (2) low-density filaments with low H2 abundance, 
and (3) high-density filaments. 

For low density filaments with xh 2 ,o ^ 3 x 10~ 3 , HD cooling over- 
whelms H2 cooling during the contraction, and the temperature de- 
scends down to T ~ 50 K accordingly. When the density reaches a 
critical density of HD (n cr ~ 10 4-5 cm -3 ), the cloud contraction tends 
to become quasistatic. Thus, the fragmentation will take place after 
that stage. For low density filaments with xh 2 ,o 10~ 3 , because 

of the low xh 2 , the temperature stays a relatively high value at a few 
hundred K, and HD cooling does not play a role during the contraction. 
The contraction proceeds quasi-statically after the density reaches the 
critical density of H2 (n > nH 2)C r = 10 3-4 cm" 3 ). As for high density 
filaments with no ^ 10 5 cm" 3 , the temperature stays at a relatively 
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Figure 1. Mass distributions of the fragment derived from ID simulations for the 
models with (a) ih 2 ,o = 1 X 10 -4 and (b) 1 x 10 -2 , respectively. The abscissa and 
ordinate denote the initial central density and the parameter /, respectively. The 
solid lines denote the contours of the fragment mass which are labeled with adjacent 
numbers. The dashed lines show the lines at which the HD cooling rate is equal to 
H2 cooling rate at the epoch of fragmentation. In the left regions of the dashed lines, 
HD cooling is more efficient than H2 cooling. 



high value because both H2 and HD cooling is saturated, and therefore 
HD does not play a significant role for the thermal evolution. The con- 
traction of such a high-density filament is accelerated by the effective 
three-body H2 formation when the density reaches 10 8-9 cm -3 . Then, 
the fragmentation is not expected to take place until the cloud becomes 
opaque to H2 lines at the density of 10 12-13 cm -3 . 

Figures la and lb show the distributions of the fragment mass de- 
rived from the ID simulations (see Nakamura and Umemura (2001b) 
for more detail). For the models with xh 2 ,o < #H 2 ,cr ~ 3 x 10 -3 (Figures 
la), the mass distribution of the fragments is quite similar to the case 
without HD (Paper II), because HD cooling does not play an important 
role in the thermal evolution of the filaments. There is a steep boundary 
at n ~ 10 4 — 10 5 cm -3 in the mass distribution of the fragments for 
/ p> 3. For the models with no ^ 10 5 cm -3 , the fragment masses take 
their minima at 1 ~ 2 M , whereas, for the models with no <^ 10 5 
cm -3 , they are greater than ~ 10 2 M Q . [For the models with low 
densities (no ^ 10 5 cm -3 ), the contraction becomes quasistatic when 
the density reaches the critical density of H2 and then linear density 
fluctuations can grow nonlinearly before the three-body H2 formation 
becomes dominant (no ;> 10 8-9 cm -3 ). In contrast, for the models 
with high densities (no ;> 10 5-6 cm -3 ), the contraction time does not 
exceed the fragmentation time until the H2 lines become optically thick 
at n ~ 10 12 — 10 13 cm -3 .] This sensitivity in the fragment mass comes 
from the rapid increase in xh 2 due to the three-body reactions. 
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On the other hand, when xh 2 ,o > 3 x 1CT 3 , HD cooling is more 
effective than H2 cooling for low-density filaments. Thus, the maximum 
mass of low-density region reduces to a few tens Mq. The minimum 
mass does not change because HD is not a dominant coolant in the 
evolution of dense filaments. Therefore, similarly to the models with 
low xh 2 ,o the dependence of the fragment mass on no exhibits a steep 
boundary around no = 10 4 — 10 5 cm~ 3 for / ^ 3. 

The existence of a steep boundary in the mass distribution implies 
that the IMF is likely to be bimodal if both low-density and high- 
density filaments bear stars. The low- mass peak is around a few Mq, 
which is not sensitive to the abundance of H2 formed in a parent cloud. 
The high-mass peak is ss 1O 2 M if xh 2 ,o ^3x 10~ 3 , while it is « 1OM 
if xh 2 ,o 3 x 10~ 3 . In the next section, we discuss some implications 
of such bimodal IMF in the very metal-deficient gas. 

4. Implications for Galaxy Formation 

Recently, Susa and Umemura (2000) investigated the pancake collapse 
of pregalactic clouds under UV background radiation. They found that 
once the pancaking disk is shielded against external UV radiation, 
the H2 molecules form efficiently via H~ reaction with abundant free 
electrons produced by UV background, and the resultant abundance 
reaches xh 2 ~ 3 x 10~ 3 . The pancake disks probably fragment into 
filaments in which stars would form. In this case, HD cooling is expected 
to become efficient in low-density filaments and then, the high mass 
peak of the IMF would go down to ~ 10 M Q . 

The time-decreasing UV background radiation is also likely to influ- 
ence star formation in dwarf galaxies which may be related to excess 
number count of faint blue objects observed in the Hubble Deep Field. 
Corbelli, Galli, and Palla (1997) studied the effects of the declined UV 
background on the thermal evolution of the protogalaxies. They found 
that there is a critical redshift of z ~ 1 — 2, below which the declined 
UV radiation is shielded by the gas disks where the H2 abundance 
reaches 10~ 2 owing to high ionization degree by the UV radiation. The 
enhanced H2 formation promotes a rapid transition toward the cold H 
I phase with ~ 10 2 K. Thus, in such dwarf galaxies, high mass peak of 
the IMF would decrease to ~ 10 M Q owing to HD cooling. 
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